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Abstract 

Using effective field tlieory we prove inequalities for the correlations of two-nucleon operators in 
low-energy symmetric nuclear matter. For physical values of operator coefficients in the effective 
Lagrangian, the S = 1, / = channel correlations must have the lowest energy and longest corre- 
lation length in the two-nucleon sector. This result is valid at nonzero density and temperature. 
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The study of inequahties in quantum chromodynamics (QCD) dates back over twenty 
years. Even though the information one may extract from these inequahties is hmited, 
they provide rare insights into the nonperturbative QCD vacuum ll|ll2lyll and hadronic 

nnnn 

spectrum [5| p [7| which currently cannot be addressed by other purely analytic methods. 
Most rigorous QCD inequalities are derived using the Euclidean functional integral. The 
close relationship to lattice QCD simulations should be apparent. Typically one compares 
the magnitude of various observable quantities by means of the Cauchy-Schwarz or Holder 
inequalities. The comparison requires positivity of the functional integral measure for at 
least one of the quantities being compared and in some cases positivity of the observable 
for arbitrary background gauge configurations. These same attributes are also important 
in lattice QCD in connection with the well-known sign/phase problem. A review of the 
literature on QCD inequalities can be found in 0|. 

Recently there has been some interest in simulations of nuclear matter on the lattice 
using effective field theory [3 [ill [l^ . These are low-energy simulations with nucleons 
and sometimes pions as point particles rather than quarks and gluons. Motivated by these 
nonperturbative studies, we consider here the subject of effective field theory inequalities 
for nuclear matter. Our focus and application will be to nuclear matter, but much of the 
interesting physics is of a universal nature and appears in other systems such as trapped 



Fermi gases near a Feshbach resonance 13| 14 1. 

It should be obvious that effective field theory does not constrain properties of the vac- 
uum, nucleon and pion masses, or other few-body spectra at zero chemical potential and 
temperature. These quantities are usually phenomenological inputs which are used to deter- 
mine the coefficients of operators in the effective Lagrangian. Hence the inequalities which 
we prove here concern correlation functions at nonzero chemical potential and/or temper- 
ature. The results are most interesting at nucleon densities where the average nucleon 
separation is less than the scattering length, and the physics is strongly coupled. 

We let represent the nucleon fields. 
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We use p{n) to represent protons(neutrons) and t(],) to represent up(down) spins. We 
use T to represent Pauli matrices acting in isospin space and a to represent Pauli matrices 



2 



acting in spin space. We assume exact isospin symmetry. In the non-relativistic limit and 
below the threshold for pion production, we can write the lowest order terms in the effective 
Lagrangian in several equivalent ways. Two of the possibilities are 

C = N[ido + £^~{m%- IJ.)]N - ICsNNNN - ICrNaN ■ NaN, (2) 



= + £;^-{m%- fi)]N - IC'sNNNN - IC'^NtN ■ NtN. (3) 

At this stage we leave out three-nucleon forces, but we will comment on this later in our 
discussion. We assume that the high-momentum modes have been removed by some regu- 
larization procedure though details are unimportant in this discussion. 

Consider states with two nucleons at the same point. In order to maintain overall 
antisymmetry, we can have either a spin S* = 0, isospin 1 = 1 state, or an S" = 1, / = state. 
These two states are eigenstates of the various interaction operators, and the corresponding 
eigenvalues are shown in Table 1. At lowest order in the pionless effective theory, we do 
not have mixing between 5" and D waves in the S = 1, 1 = channel. 



Table 1: Eigenvalues of interaction operators 





S = 1,1 = 


5 = 0, / = 1 


i : NNNN : 


1 


1 


i : NaN ■ NaN : 


1 


-3 


i : NtN ■ NtN : 


-3 


1 



(4) 



Clearly 

: NNNN := -| : NaN ■ NaN : -i : NtN ■ NtN :, (5) 



C'u — —Ct, C'g — Cs — 2Ct- (6) 

In the real world two-nucleon low-energy scattering is attractive in both the S = 1, 1 = 
and S = 0, 1 = 1 channels. It is also more strongly attractive in the S = 1, 1 = channel. 
From this we deduce that 

Cs < SCt, Ct < 0, (7) 
C's<-Ci, Ci>0. (8) 
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These inequalities can be verified nonperturbatively by solving the Schrodinger equation 



and determining phase shifts 



111 ] or by summing all bubble diagrams 21[. The analysis is 



straightforward since we have only a simple contact potential. 
The grand canonical partition function is given by 

Zg oc [ DNDNexp{-SE) = [ DNDNexp i I d^xCs 



(9) 



where we choose to write Ce 



Ce = -N[d4 



2m M 



+ {m% - ij)]N - ^CgNNNN - \CuNfN ■ NfN. 



(10) 



Using Hubbard-Stratonovich transformations [22|[2J|, we can rewrite Zq as 



Zgoc / DNDNDfD(l)exp[ I (fxC^i" 



(11) 



where 



= -N[d, -^ + (m%- i,)]N + C'sfNN + iC'J- NfN - V(/, $) (12) 



and 



-V(/, 



(13) 



We note that < and > 0, and so the / and integrations are convergent. 

Let M be the nucleon matrix generated by the background Hubbard-Stratonovich fields, 



M = - 



We note that 



T2MT2 = M* 



(14) 



(15) 



where M* is the complex conjugate of M, not the Hermitian adjoint. We conclude that 
det M must be real. Since M does not couple up and down spins we can write M with two 
identical diagonal blocks, 

M| 



M 



M 



Mt = M|. 



(16) 



We note that det M| must also be real, and so 



det M = det det = (det M^f > 0. 



(17) 



4 



Although not relevant to the main line of this discussion, we point out that actually det > 
0. This is because T2 is antisymmetric and therefore the real eigenvalues of M| are doubly 
degenerate 

We note the similarity between (jl5j) and the relation 

75M75 = Mt (18) 

for the quark matrix in QCD at zero chemical potential. From (jl8j) one can show positivity 
of the quark matrix determinant for two degenerate flavors and the dominance of pion 
correlation functions [2^. In our case, however, we observe that (fT3j) is also valid at 
nonzero chemical potential. 

Let us now consider the two-nucleon operator 

A{x) = [NMr2h[NiUx) (19) 

where N-^ and Ni are up and down spin projections of A^. Since A has quantum number 
/ = 0, overall antisymmetry requires it also have quantum number S = 1. The two-point 
correlation function for A is 

(A(x)At(0)>^^^ = {[NMr,],,[N^],ix) [NlUr,U[N;UO))^^^. (20) 

Using our Euclidean functional integral representation, we have 

{A{x)A\0))^^^ = J De[M^\x,0)UT2]t4M^\x,0)U[r2hj, (21) 

where DQ is the positive normalized measure defined by 



DfDcl) det M exp - / d^x V(/, <; 

DQ = \ (22) 

/ DfD$ det M exp ( - / d^x V(/, ^ 



From the symmetry properties of T2 under transposition and (fT3j) . 

[M^\X, 0)]jl[T2]lm[M^'\x, 0)]im[T2]ij = [M^\x, 0)],, [Ts],^ [Mf ^ (x, 0)]im[T2]ml 

= [M^\x,0)UM^~\x,0)Y,i 

|2 



Thus 



5^|[M-i(x,0)],.f . (23) 



{Aix)A\0))^^^ = I De5^|[M-i(x,0)],.|'. (24) 
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Consider now the more general operator 

B{x) = [NMbUNiU^) (25) 
for arbitrary 2x2 matrix [b]ij. We find 

{Bix)B\0))^^^ = I DQGb{x), (26) 

where 

Gb{x) = [M;\x,0)]Abrrni[M{\x,0)U[b],,. (27) 
From the Cauchy-Schwarz inequahty, 



ijlm y i'j'l'm' 
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m'l' 



<5^|[Mfi(x,0)].^|'E|[%.r (28) 

im jk 

We conclude that 

B{x)B\0))^^^\ < I De\GB{x)\<j{A{x)A\0))^^^ (29) 
where 

l = J2\i^hkf. (30) 

jk 

Let us take the zero temperature limit at fixed nonzero nucleon density. This will require 
tuning the chemical potential ^ to match the desired density as we lower the temperature. 
If we take the separation x ^ oo in the Euclidean time direction, we deduce from (j^^ that 
the energy of the lowest state with the quantum numbers of A (i.e., 5* = 1, / = 0) must be 
less than or equal to the lowest state with the quantum numbers of B, 

El<El. (31) 

This result was known from the outset for zero density. The non-trivial result is that it 
also holds true at nonzero density. In particular if we are in a region of the phase diagram 
where the deuteron is unbound, then all states with the quantum numbers of B must also 



be unbound. A formula due to Liischer 



tells us that in a periodic box of length L, the 



energy of the lowest two-particle scattering state relative to threshold is 



where is a is the scattering length with sign convention, 

1 1 

fccot5o = — + -rofc2 + .... (33) 

By combining ()31|) and ()32|). we find in the unbound case the following inequality for the 
scattering lengths Q], 

aA<aB- (34) 

If we take the limit x — > oo in any spatial direction, we find that the correlation length 
for A must be greater than or equal to the correlation length for 5, 

U > is- (35) 

This holds true for any density and any temperature. In particular if diverges, as would 
occur in a superfiuid phase transition, then must also diverge. In that case one can use 
(|29|) and (j30|l to bound the relative sizes of the dinucleon condensates. 

The inequalities ()31|1 . (jBl}, and were proven for operators B of the specific form given 
in ()25j) . It is straightforward to generalize to the case when the nucleons have arbitrary 
spin directions and/or are not at the same point. Let 

B,^y{x) = [N{x + y)m,,[N{x - y)],, (36) 

where \b]ij is a general 4x4 matrix. One can show for some positive parameters 71, 72, 73, 

74, 



lj.,T 



< 71 {A{x - y')A^{-y)) + 72 {A{x - y')A\y)) 

+ 73 {Aix + y')A\-y)) + 74 (Aix + y')A\y)) . (37) 



From ()37p we conclude that the inequalities ()31|) . ()34|) . and ()35|) hold for any two-nucleon 
operator B, 

B{x) = [ d'y [N{x + yMh{y)UN{x - y)]„ (38) 
Jn 

where the region of integration Vt is bounded. 

In summary we have shown that the = 1, / = channel must have the lowest energy 
and longest correlation length in the two-nucleon sector. These results are valid at arbi- 
trary density and temperature and can be readily checked in lattice simulations using chiral 



effective field theory. Assuming that the effective Lagrangian can describe the relevant 
physics, it would be interesting to study two-nucleon correlations near the region of the 
phase diagram where the deuteron becomes unbound. In principle these results can also be 
checked in lattice QCD at finite density, though the complex phase problem will make the 
calculation exceedingly difficult. 

It turns out that the addition of the three-nucleon force term {NNY j^l can be intro- 
duced by modifying the coupling of / and NN while preserving positivity of det M and 
the conjugate similarity relation in ()15|). Other modifications can also be made for various 
higher order terms in the effective Lagrangian. These and other extensions will be discussed 
in future work. 

The author thanks Paulo Bedaque, Jiunn-Wei Chen, and Thomas Schaefer for several 
helpful discussions. He also thanks Paulo for organizing the Summer of Lattice Workshop 
2004 at Lawrence Berkeley Laboratory where part of this research was completed. This 
work was supported by the Department of Energy. 
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